Fano effect of strongly interacting quantum dot in contact with superconductor 
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The physics of a system consisting of an Aharonov Bohm (AB) interferometer with a single level 
interacting quantum dot (QD) on one of its arms, and attached to normal (N) and superconducting 
(S) leads is elucidated. Here the focus is directed mainly on N-AB-S junctions but the theory is 
capable of studying S-AB-S junctions as well. The interesting physics emerges under the conditions 
that both the Kondo effect in the QD and the the Fano effect are equally important. As the Fano 
effect becomes more dominant, the conductance of the junction is suppressed. 



Motivation: Transport through a mesoscopic 
Aharonov-Bohm (AB) ring (with an interacting 
quantum dot (QD) situated on one of its arms) weakly 
attached to normal (N) metallic leads is the subject of 
recent intensive experimental Q, 1^ and theoretical 
HIE IE HIE studies. To the above class of experi- 
ments one may also add STM measurements on a single 
magnetic atom adsorbed on a gold surface ^E|- The 
main result of the theoretical analysis |E indicates 
that in such N-AB-N junctions there is an interplay 
between two fundamental physical phenomena, namely, 
the Fano and the Kondo effects. The Fano effect is 
related to an interference between two electron waves one 
passing through the quantum dot (with a discrete level) 
and one travelling along the direct 'reference' channel 
characterized by its continuous spectrum. It results in 
an asymmetric shape of the conductance as function of 
the applied bias (or gate) voltage. The Kondo effect is 
one of the simplest manifestations of many-body physics 
exhibiting strong correlations. It plays an important 
role in electron transport through quantum dots, where 
the role of a magnetic impurit y is played here by the 
presence of localized electrons [ll|, |l^- The interplay 
between these two effects causes the suppression of 
Kondo plateau with increasing transmission through the 
direct channel. 

An analysis of the conductance of an AB interferom- 
eter when one of the leads is a normal metal and the 
other one is a superconductor (S) is the subject of the 
present research. The geometry of such N-AB-S junction 
is schematically displayed in FigC] We derive a general 
formula for the conductance of N-AB-S and S-AB-S junc- 
tions. 

As for the motivation, there is indeed a grow- 
ing interest in these types of superconducting junctions, 
especially in the light of recent experiments which in- 
volve carbon nanotubes as weak links |l^ Hsf. They 
are characterized by relatively high Kondo temperature 
which helps the formation of Kondo plateau in conduc- 
tance experiments. The Kondo finger-prints in such junc- 
tions is expected to survive even for cases where the su- 
perconducting gap approaches the Kondo temperature 
from below, that is, A - (see Refs.|lE|ll|l3)- 



FIG. 1: Schematic form of an N-AB-S junction. 



Model Hamiltonian and the current: In the system 
under consideration (see Fig. transport of an elec- 
tron between the two-dimensional electrodes (N on the 
left and S on the right) is possible via two paths, ei- 
ther through the QD, or else through the other 'direct' 
channel. The dynamics of the system is governed by the 
Hamiltonian 



H = Hl + Hb + Hd + Ht -\- Hlr, 



(1) 



in which iJj, (j = L,R) are the Hamiltonians of the 
electrons in the electrodes which depend on the electron 
field operators V'jcr(r,t) where r = {x,y) and cr = ± is 
the spin index. 



H, 



dr[V'],(r)e(V)^,. - aV^]^(r)V']^V',i(r)^,T W].(2) 



Here a is the BCS coupling constant and ^(V) = 
— V^/2m — /X with /i being the chemical potential at tem- 
perature T. The dot is represented by a single level An- 
derson impurity with energy eo < and Hubbard repul- 
sion U. Later on we will concentrate on the Kondo regime 
by setting J7 — > oo and assuming that |eo| exceeds any 
other energy scale except U. The Nambu representation 
for all fermion operators {ipa- (r) for the lead electrons and 
Ca for the QD electrons) is employed below (for example 
= (c|,C2) — (ct,C|)). The dot Hamiltonian then reads. 



(3) 
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where n^- = cl-Ca and the Pauh matrix Tz acts in Nambu 
space. Assuming a symmetric junction, the tunnelUng 
Hamiltonian Hf from the leads to the dot is, 



' h.c, 



(4) 



where T is the tunneUing amplitude. The occurrence of 
direct transport channel is represented by the term Hlr 
between left and right leads. The gauge is chosen such 
that the AB flux appears only in this direct term, 



LR 



WV^t^(0)T,FV'i(0) + /i.c., 



(5) 



where F = exp(z</3T2) is the AB phase factor. Note that 
unlike the case of N-AB-N junctions, the phase factor is 
non-Ahelian. 

In the mean field slave boson approximation (MFA) 
which is employed below to handle the Kondo problem 
(in the limit C/ — + oo), the tunnelling amplitude is mod- 
ified by replacing the (slave) boson operators by their 
mean-field values, (6 ^< h >) and the Hubbard repul- 
sion term is dropped. This procedure imposes a con- 
straint on the number of bosons and fermions and the 
Hamiltonian must now include a term which prevents 
double occupancy in the limit U ^ oo, 



(6) 



where is the renormalized new level position which has 
the role of a Lagrange multiplier. 

Formula for the current at one lead, say the right 
(supreconducting), is derived from the basic relation 
I = —ie < [nn, H] >. It is convenient to display the 
result of the commutation relation in terms of Keldysh 
Green's functions (GF) , namely. 
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\^Tr{F^G'tj, + H.c) + ^TpTriG'^j^ + H.c)\7) 



where p is the density of states in the leads (assumed 
to be the same in both), 7 = Tr^W^p^ and the trace 
operation acts in Nambu and energy spaces (Tr[A(a;)] = 
/ Aii{uj)(Lj). Other notations are Gf^ = Glr~^^lr 
(similarly for G^j^), where G^j^{t) = i < ^/'JjV'L(i) > and 

G^nit) = ^ < Vlic(t) >. 

Employing equation of motion method for all Keldysh 
GF it is possible to express the current in terms of the dot 
GF. We single out the current through the direct channel 
by writing Gf^ = G{ 



'LdR 



with. 



= 2^[gLLFTzgnr (8) 
G^R - n9LLil + 2V^FTzgn)Gil + 2^gLFTz)gRj,f 
Gf« = T[TzGil + 2^gLFrz)gRRf. (9) 

Here F = 27rpT^ stands for tunneling rate through the 
dot, G = TzGtz, G is the dot GF, gL,R are generally 
nonequilibrium GF of the leads (left, right). The other 



two GFs include multiple scattering events and have the 
form, 



gLL (1 - ■^ighFTzgRTzF^) ^gL 
9RR = (1 - '^IQRF'^TzgLTzFy^gR, 



(10) 



Each GF {gL,R, gLL, Qrr, and G) has a standard 2x2 
Keldysh matrix structure with elements Gn = G^ , 
G22 = G-^, G21 = 0, G12 = G^ . The superscript K 
stands for Keldysh component of the matrix product oc- 
curring in the square brackets of Eq.(7-|H). The set of 
equations (T-^J together with the expression (0 for the 
current formally complete our task . They give an expres- 
sion for the current through a closed AB interferometer 
in terms of the GF of the (strongly interacting) QD. The 
above formalism has been tested for the case of linear 
conductance with normal leads on both sides (N-AB-N 
junction) and the result agrees with pertinent calcula- 
tions 

These rather general expressions are now employed for 
elucidating a particular case of interest, namely, the zero 
temperature limit of the linear conductance in N-AB-S 
junctions. Specifically, as in Fig. the left lead is a 
normal metal biased with an external voltage V whose 
GFs are, ^ = - ± i/2, g^ = -i(tanh(w/2T) + S/tz); 
Sf = ey/2T(cosh(w/2T))-2. On the other hand, the 
right lead is an unbiased (s-wave) superconductor with 
gap A whose equilibrium GFs have a standard formjlSl]. 

Mean Field Approximation and Conductance: Within 
the MFA the algorithm starts with calculations pertain- 
ing to a geometrically identical system albeit with nonin- 
teracting QD (J7 — > 0). Then, at the end, the replacement 
r — > , e — + is executed in the conductance formula. 
These two quantities (the Kondo temperature Tk = 
and the effective renormalized position of the level) are 
evaluated by solving two mean field equations similar to 
those derived in Refs. At zero bias one is free to 

adopt the Matsubara form of these equations. The first 
one, 27r(eT. — e)-|-rrr(t/A) = 0, follows from the extremum 
requirement of the effective action over Cr- The second 
one, 2-kTk + TTr{QTz) — 0, reflects the single occupation 
condition. Here 

A(w) ^iuj- ErTz + {tz ~ iy/^ F)g rii{Tz - iy/ly^F'') 

\UJ\ \UJ\ 



Q ^{uj) = iuj - erTz - TkA, 



(11) 



^--^^^ = -2dR(- + 7^7^:^- .r.A), (12) 



where d{Lu) = (1 + 7^))w^VA2 + u'^ + 2jll>^\uj\. Expres- 
sions (|11I12|) for the dot Matsubara GF were obtained 
from the effective action derived from the Hamiltonian 
of the N-AB-S system. The self-consistence equations 
were solved for different values of the direct transmission 
parameter 7 while the level positions e are scanned over 
the wide range (-0.4,0.1). The parameters are tuned such 
that the carbon nanotubes in Kondo regime considered 
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in reference [13 can be satisfactorily regarded as an An- 
derson impurity. The half width of the normal electrode 
density of states is served as an energy unit and the tun- 
nelling rate through the dot is fixed at F = 0.15. The 
superconducting gap A is set at a value A = 0.03. For 
the majority of level positions e the A < Tr- and thus 
the MFA is hence justified. 

To obtain the zero bias differential conductance for the 
noninteracting dot, expressions for the GF G^,G^,G^ 
are found for this case and are then directly substituted 
into equations H7I9|) . It is important to note that any 
element in these GF which is related to the superconduc- 
tor electrode must be evaluated at energies |u;| < A (this 
is not the case for the self consistence equations). For 
a dot modelled by a simple resonance energy level, the 
formulae for the GF is (compare Eg. 1)11(1 '). 

G = W - er, -T[{t, + 2^gLF)gRRiT, +2^gLF^)]}-' 

Its retarded component reads then, 

= [9 + Tz9z+Txgx + Tygy]~'^, (13) 

where g — ui + r|(l — P(l — 7)), the coefficients of 
the Pauli matrices are = — e — P cos{ip) , g^ — 
-irQ(l 7COs(2(/?)) and gy = irQ7 sin(2.^). On the 
Fermi surface the relations P — —7/(1 + 7'^), Q = 
1/(1 -1-7^) hold. Inserting these GF into equations H7I9|) 
for the current one arrives at an expression for the zero 
bias conductance. After lengthy calculations the conduc- 
tance is obtained and checked to be an even function of 
the flux, as is dictated by Onsager relations (although it 
is not immediately apparent, as asymmetric terms can- 
cel). 

2 3—^ 

aNs = ^{Tw+'^mv)Y.^jCos{j^)) (14) 

3=0 

Here Tw = 47/(1 + 7^)^ is the background (direct) NS 
transmission. Introducing the quantity T/vat = 47/(1-1- 
7)^ which is the analogous quantity (direct transmission) 
for the N-AB-N junction^ then Tw = / {2 - Tn n)'^ ■ 
This is precisely the relation between transmission coef- 
ficients for N-BB-N and N-BB-S junctions (here BB is a 
"black-box" representing any non-interacting scatterer) 
suggested in Ref. [l9i| . using Landauer scattering matrix 
approach. The coefficients T\y and Tj in equation l|14|) 
do not depend on flux. Note, however, that the normal- 
ization factor M is an even function of the flux, 

[AA((^)]-i/2 = 2(l + 72)3[r2(l + 7) + 2e2(l+72) 
-4re73/2 cos{(p) + r'^'y cos(2(^)] 

To = -8r[r2(-i + 272 + 87=^ + 57*) + 

167e2(-l + 7-72 + 73 + 274)] 
Ti = 64V7e[r^(-l-27^ + 

T2 = -64r7[rV(l+7) + e^(-l + 57*)] 



The last two coefflcients are T3 = 128^/77'^r^e and 
T4 = — 167^r^. Here the result of calculations of con- 
ductance as function of gate voltage p4|l (with F — 0.04 
and zero AB phase) is shown on Fig[21 A typical Fano 
asymmetry form in the case 7 = 0.1 is easily detected. 
As was noted above, in the MFA, F and e appearing in 




FIG. 2: Linear conductance of an N-AB-S junction (see Fig. 
[g at r = and zero magnetic field as function of level posi- 
tion e at 7 = 0.1 (solid line) and 7 — 0(dash line) 

the above equations should respectively be replaced by 
Tk and e^, which, in turn, are obtained through the solu- 
tions of the self-consistence equations. The conductance 
of an N-AB-S junction as function of the level position 
e is displayed in flgure |2| and as a function of the flux in 
flgure 01 



h 




FIG. 3: Linear conductance of an N-AB-S junction (see Fig. 
dl at r = and zero magnetic field as function of level posi- 
tion e at 7 = 0; 0.1 and 7 = 0.2. 

Discussion: If direct tunnelling is completely suppressed 
(7 — 0) the electron trajectory passes solely through the 
quantum dot. This is a single channel Coulomb block- 
ade situation for an N-QD-S junction. For level energies 
which satisfy pure Kondo hmit, |e| >> F, one expects 
Kondo behavior and Kondo plateau in the differential 
conductance at zero bias^3- We performed calculations 
for this range of e and found such a dependence of the 
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This property characterizes solely the Kondo physics. It 
has also been proved by performing the calculations in 
the strict Kondo limit |e| >> F. For the higher energy 
in Fig. 4, the QD is in the mixed valence regime and the 
conductance maximum is not pinned at ip — 7r/2. 

This research is partly supported by grants from the 
Israeli Science Foundation (ISF) and the American Israeli 
Binational Science foundation (BSF). We would like to 
thank Moshe Shechter and B. Bulka for very helpful dis- 
cussions. 



FIG. 4: Linear conductance of an N-AB-S junction (see Fig. 
at r = versus AB phase for different level positions {e—- 
0.4 top line, -0.34 middle line, -0.24 bottom (red) line. All 
curves correspond to 7 = O.f . The parameters are tuned to 
cover the Kondo and crossover to Kondo regimes. 



linear conductance. If the above inequality is less pro- 
nounced (as is typical for carbon nanotubes) , the Kondo 
effect is not completely destroyed: The linear conduc- 
tance approaches the unitary limit as |e| grows. In fig- 
ure 3, depicting the conductance as function of the level 
position at zero AB flux and different background trans- 
missions (T/vAT = 0,0.33,0.55, or 7 = 0,0.1,0.2) clearly 
shows this behavior. In case 7^0 there is a clear sup- 
pression of the Kondo effect and the attenuation of the 
conductance at its plateau step is evident although, in 
this region of energies, the Kondo effect survives. This 
also follows from solution of the mean field equations (for 
7 = 0.1) yielding an effective Kondo temperature Tk 
which is somewhat smaller than in the case 7 = 0. The 
situation here is reminiscent of the Fano Kondo effect in 
normal junctions 0. Such behavior of the conduc- 
tance may serve as an indication of Kondo correlations 
in the dot. 

The crossover to the Kondo regime reveals itself more 
clearly in Fig. 4. It manifests a remarkable behavior of 
the conductance as function of the magnetic field (the AB 
flux) for three different values of the dot level position. 
For the two lower energies there is a maximum very close 
to the Unitary limit 4e^//i which is pinned near ip — tt/2. 



[1] A.Yacoby, M. Heiblum, D. Mahalu, and H. Shtricman, 

Phys. Rev. Lett. 74, 4047 (1995); Y. Ji,et al. Science 

290, 779 (2000). 
[2] W. G. van der Wiel et al, Science 289, 2105 (2000). 
[3] K. Kobayashi, H. Aikawa, S. Katsumoto and Y. lye, 

Phys. Rev. Lett. 88 256806 (2002). 
[4] U. Gerland et al. , Phys. Rev. Lett. 84, 3710 (2000). 
[5] J. Konig and Y. Gefen , Phys. Rev. Lett. 86, 3855 (2001). 
[6] W. Hofstetter, J. Konig, and H. Schoeller , Phys. Rev. 

Lett. 87, 156803 (2001). 
[7] B. R. Bulka and P. Stefanski, Phys. Rev. Lett. 86, 5128 

(2001) . 

[8] M. A. Davidovich et al, Phys. Rev. B 55, R7335 (1997). 
[9] O. Entin-Wohlman et al, Phys. Rev. Lett. 88, 166801 

(2002) . 

[10] V. Madhavan et al Science 240, 567 (1998). 

[11] L. I. Glazman and M. E. Raikh, Pis'ma Zh. Eksp. Teor. 

Fiz. 47, 378 (1988) [JETP Lett. 47 452 (1988)]; T. K. 

Ng and P. A. Lee, Phys. Rev. Lett. 61, 1768 (1988). 
[12] Y. Meir, N. S. Wingreen, and P. A. Lee, Phys. Rev. Lett. 

70, 2601 (1993). 
[13] M.R. Buitelaar , N. Nussbaumer, and C. Schonenberger 

Phys. Rev. Lett. 89, 256801 (2002). 
[14] A. Yu. Kasumov et al. Science 284, 1508 (1999). 
[15] J. Nygard, D. H. Cobben, and D. E. Lindelof, Nature 

408, 342 (2000). 
[16] A. A. Clerk, V. Ambegaokar, and S. Hershfield, Phys. 

Rev. B 61, 3555 (2000). 
[17] P. Schwab and R. Raimondi, Phys. Rev. B 59, 1637 

(1999); R. Fazio and R. Raimondi, Phys. Rev. Lett. 80, 

2913 (1998); 82, 4950 (E)(1999). 
[18] Y. Avishai, A. Golub, and A. D. Zaikin, 

|cond-mat/0111442| Phys. Rev. B67(RC), 041301 

(2003) . 

[19] C. W. J. Beenakker Phys. Rev. B 43, 134515 (1991). 



